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1 Introduction 



At the beginning of 2011 the status of the inclusive cross-section for Higgs-boson production in gluon 
fusion was summarized [1 . Corrections arising from higher-order QCD, electroweak effects, as well as 
contributions beyond the commonly-used effective theory approximation were analyzed. Uncertainties arising 
from missing terms in the perturbative expansion, as well as imprecise knowledge of parton distribution 
functions, were estimated to range from approximately 15—20%, with the precise value depending on the 
Higgs boson mass. For an updated study we refer to Ref. [2J. 

Recently the problem of going beyond the zero-width approximation has received new boost from the 
work of Refs. [314] and of Ref. [5] which implemented the complex-pole-scheme with an estimate of the 
residual theoretical uncertainty (see also the work of Ref. [5]). Here we only recall that the complex pole 
describing an unstable particle is conventionally parametrized as 



with i = W, Z,H etc. (see Ref. [51715] 1. 

In the current experimental analysis there are additional sources of uncertainty, e.g. background and 
Higgs interference effects [9110111112113] . As a matter of fact, this interference is partly available and should 
not be included as a theoretical uncertainty; for a discussion and results we refer to Refs. [14115116117] . In 
particular, from Refs. [T3] we see that at nu — 600 GeV (the highest value reported) the effect is about 
+40% in the window £ = 440—560 GeV, where £ is the Higgs virtuality; the effect is practically zero at the 
peak and reaches —50% after C, = 680 GeV (no cuts applied). For the total cross-section in gg — > lvlV at 
/ih = 600 GeV the effect of including the interference is already +34% and rapidly increasing with fin. 

We stress that setting limits without including the effects of the interference induces large variations in 
rate and shape that will propagate through to all distributions. Therefore, any attempt to analyze kinematic 
distributions which are far from the Standard Model (SM) shape may result in misleading limits. 

The importance of a complete understanding of the shape of the background from non-resonant diagrams 
has been emphasised in Refs. 18 19 . It was shown in this work that a heavy Higgs boson with mass larger 
than 800 GeV does not lead to a pronounced peak structure in the lineshape and predictions for the non- 
resonant background must therefore be as accurate as possible in order to discriminate a heavy Higgs boson 
from a light one. 

In the current experimental analysis for heavy Higgs searches, a theoretical uncertainty of 150 /^[%] (/in 
in TeV) has been used for conservative estimat^. For a Higgs boson of 700 GeV this amounts to ±51%. 

One might wonder why considering a Standard Model (SM) Higgs boson in such a high-mass range. 
There are classic constraints on the Higgs boson mass coming from unitarity, triviality and vacuum stability, 
precision electroweak data and absence of fine-tuning [20J . The situation is different if we consider extensions 
of the SM: in the THD model, even if the SM-like Higgs boson is found to be light (< 140 GeV), there is 
a possible range of mass splitting in the heavy Higgs boson. In general, for a given Higgs boson mass, 
the magnitude of the mass splittings among different heavy scalar bosons can be determined to satisfy the 
electroweak precision data, see Ref. [21] . 

This paper is organized as follows. In Section[2]we discuss results available in the literature. In Section[3] 
we present and discuss numerical results. Conclusions are presented in Section 3] 

2 Summary of available results 

There are serious problems in including the signal/background interference in gluon-gluon fusion and 
very few examples of theoretical predictions, e.g. interference has been computed for the di-photon signal 
in Ref. [22]. Let us concentrate on the process gg — > TIL, the whole cross-section can be written as follows 
(here, for simplicity, we neglect folding the partonic process with parton distribution functions): 



Mi -ifM7%, 



(1) 



Cgg^zz - (T gg ^. zz (S) + er gg ^ zz (/) + crgg^zz(-B), 
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where S, B and I stand for signal (gg — > H — > ZZ), background (gg — ► ZZ, i.e. gluon initiated box 
contribution for all three quark doublets) and interference; the signal can be written as 

1 C 2 r H ^zz 



' : : -Zz(S) — O-gg^H^Zz(C) — — CT gg^H 2 77 ' ( 3 ) 

C - «H 
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where we have introduced sh, the Higgs boson complex pole |7I5] ; furthermore, £ is the Higgs boson 
virtuality. So far, most if not all theoretical prediction have been devoted to compute the signal with the 
highest possible precision, next-to- leading order (NLO) and beyond. Therefore, in Eq.@, we have 

• the production cross-section, <7 g g_>H , with next-to-next-to- leading logarithmic resummation (NNLL) |23j 
i.e. with NNLL + NNLO + EW + bottom quark contribution up to NLO + NLL (three loop level 
plus resummation), see Ref. [5] 

• the partial decay width of an off-shell Higgs boson of virtuality £, Lh^zz, which is known at NLO 
(one-loop) with leading NNLO effects in the limit of large Higgs boson mass, see Ref. |24l25j . 

However the background (continuum gg — \ ZZ) and the interference are only known at leading order (LO, 
one- loop) [26 . Here we face two problems, a missing NLO calculation of the background (two- loop) and the 
NLO or NNLO signal at the amplitude level, without which there is no way to improve upon the present LO 
calculation. For low values of the Higgs boson mass the interference arises primarily from the imaginary part 
of the continuum background interfering with the real part of the signal. For fin < 2 mt , where we can use 
the effective theory (i.e. large- mt limit), it would be relatively easy to get the signal amplitude; however, 
what might be the tougher part is implementing it in the same program that contains the background. 
Above the tt -threshold, getting the signal amplitude becomes more difficult. We know that the effective 
theory misses imaginary parts in this region, and it is not clear how one would trust a calculation for the 
interference using it. 

Of course, putting in the O (a™) corrections to the signal without the background can only be considered 
an approximation to the interference correction. It is difficult to believe that background NLO calculation 
will be done in a foreseeable future since not all basic master integral for the two-loop contribution are 
known, at least analytically. 

We can also summarize the basic features of the LO interference, especially for high value of the Higgs 
boson mass where the contribution from the imaginary part of the signal is not negligible. The main issue 
is on unitarity cancellations at high energy. Of course, the behavior of both LO amplitudes (signal and 
background) for Mzz — > oo is known and simple and any correct treatment of perturbation theory (no 
mixing of different orders) will respect the unitarity cancellations. Since the Higgs boson decays almost 
completely into longitudinal Zs, for Afzz — > oo we have |26j (for a single quark q) 

M2 7 m 2 n „ M|, m 2 n , M|„ 

^s-^F^AHln 2 ^ A B ~- — ^ ln 2 ^, 4 
2M^ m^ 2M£ m z q 

where Ah is the Higgs propagator, showing cancellation in the limit (M| z Ah — >• 1). However, the behavior 
for M| z — > oo (unitarity) should not/cannot be used to simulate the interference for Mzz < Mh- The only 
relevant message to be derived here is that unitarity requires the interference to be destructive at large values 
of Mzz- The explicit LO calculation also shows that the interference is constructive below the Higgs peak. 
The higher-order correction in gluon-gluon fusion 27 28 29 30 31 32 33 23 have shown a huge K -factor (for 
updated cross-sextions at 8 TeV see Ref. [34 ) 

NNLO 

K = ' CT Pr°d = ^gg^H- (5) 

^"prod 

A potential worry, already addressed in Ref. [141 . is: should we simply use the full LO calculation or should 
we try to effectively include the large (factor two) if -factor to have effective NNLO observables? There are 
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different opinions since interference effects may be as large or larger than NNLO corrections to the signal. 
Therefore, it is important to quantify both effects. So far, two options have been introduced to account for 
the increase in the signal (c/.Ref. [3S 36 ). Let us consider any distribution D, i.e. 

da da 
or 

dMzz dp\ 

where M zz is the invariant mass of the ZZ -pair and p\ is the transverse momentum. Two possible options 



D =^7T- or 3Tz etc - ( 6 ) 



additive where one computes 
multiplicative where one computes 



^NNLO _ ^NNLO (5) + £)LO (J) + £,LO (_g) ^ 



) e T° = K D [D™(S) + D LO (I)] + D-°(B), K D = ^ of { > , (8) 



where K rj is the differential K -factor for the distribution. 



In both cases the NNLO corrections include the NLO electroweak part, for production [37] and decay [25] . 
It is worth noting that the differential if -factor for the ZZ -invariant mass distribution is a slowly increasing 
function of M zz , going (e.g. for = 700 GeV) from 2.04 at M zz = 210 GeV to 2.52 at M zz = 1 TeV. 

The two options, as well as intermediate ones, suffer from an obvious problem: they are spoiling the 
unitarity cancellation between signal and background for M zz — > oo, breakdown which is described in 
details in Ref. |26j . Therefore, our partial conclusion is that any option showing an early onset of unitarity 
violation should not be used for too high values of the ZZ -invariant mass. 

Therefore, our first prescription in proposing an effective higher-order interference will be to limit the 
risk of overestimation of the signal by applying the recipe only in some restricted interval of the ZZ -invariant 
mass, e.g. [/in — 7h , Mh + 7h]- This is especially true for high values of /in where the width is large. 

Explicit calculations show that the multiplicative option is better suited for regions with destructive 
interference while the additive option can be used in regions where the effect of the interference is positive, 
i.e. we still miss higher orders from the background amplitude but do not spoil cancellations between signal 
and background. 

Actually, there is an intermediate options that is based on the following observation: higher-order cor- 
rections to the signal are made of several terms (see Ref. [29] for a definition of Act) , 



4 , J i 4 Jzq Jz U 



hi h3 



where the sum is over incident partons; furthermore £ = z s is the Higgs virtuality, zq is a lower bound on 
the invariant mass of the H decay products, K = v s is the invariant mass of the incoming partons and 
the luminosity is defined by 

1,1 dx 



I dx / u \ 

£ij(v)= — fi(x,fi F ) fj (-,/*) ■ (io) 

./ „, X V X ' 



The partonic cross-section is defined by 

>H+x(C> K , MR, Mf) = CTgg^H 5 (l - - J + - (A(T gg ^Hg + Acr qg ^Hq + A(Tq q ^Hg + NNLO) . (11) 



E 



From this point of view it seems more convenient to define 



, KS £, rarest ^gg D NNL ° (gg ~> H(g) -> ZZ(g)) 

K» = K U +K D , d = ^o (gg ^ H ^ZZ) (12) 



and to introduce a third option 
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• intermediate 

^nnlo = Ru d lo^ + ( K ^) 1/2 d*°(I) + D W {B) (13) 

which, in our opinion, better simulates the inclusion of K -factors at the level of amplitudes (although 
we are still missing corrections to the continuum amplitude). 

Alternatively one could consider a different approach when M zz > > /i H ; it is based on Eq. and on the 
work of Ref. [35] and amounts to neglect the background (where NLO corrections are not available) while 
modifying the Higgs propagator, 

1 =(l + i l!L)(M* z -mI+ 1 E^a4 z Y 1 -> ^(Mlz-Ml+il^-MizY 1 (14) 

M| z - s H V M H A zz h M h ' M$ z I zz M H / 

where we have introduced mass and width in the Bar-scheme jS] according to 

^I=Mh+7h ^Hr H =A/H7H. (15) 

This recipe, Eq. (fT4|) . should be used only for M zz >> /in and should not be extended below the resonant 
peak. 

In the following Section we present numerical results in the high Higgs-mass region. 

3 Numerical results 

In the following we will present numerical results obtained with the program HTO (G. Passarino, unpub- 
lished) that allows for the study of the Higgs-boson-lineshape, in gluon-gluon fusion (ggF), using complex 
poles. HTO is a FORTRAN 95 program that contains a translation of the subroutine HIGGSNNLO written by 
M. Grazzini for computing the total (on-shell) cross-section for Higgs-boson production (in ggF) at NLO and 
NNLO 39 40 41] and a translation of the program ggzz by E.W Glover and J.J. van der Bij for computing 
the LO interference in gg — > ZZ. 

All results in this paper refer to ^/s = 8 TeV and are based on the MSTW2008 PDF sets g3]. They are 
implemented according to the OFFP - scheme, see Eq.(45) of Ref. JSJ . Furthermore we use renormalization 
and factorization QCD scales that evolve with the Higgs virtuality (Mzz). 

In Table[T]we present the effect of the interference w.r.t. signal -I- background for the total cross-section. 
We select a leptonic final state (the branching ratio for both Z bosons to decay into e or [i is 4.36 10 -3 ) and use 
p z > 0.25 Mzz and 2Mz < Mzz < 1 TeV. As is evident the strong cancellations between the constructive 
and destructive interference below and above the peak result in a small effect on the total cross-section. 
However, the effect is drastically different on distributions. We present results for the Mzz -distribution 

Table 1: Interference effect for the total cross-section at LO and NNLO with multiplicative (M) option of Eq.([8]), 
additive (A) option of Eq.0 and intermediate option of Eq.|T3]). Here p| > 0.25 Mzz and 2 M z < Mzz < 1 TeV. 
In the table we show the percentage effect of the interference, i.e. interference/ (signal + background). 

/i H [GcV] LO NNLO(A) NNLO(I) NNLO(M) 

400 ' 0.80[%] 0.64[%] 1.05[%] 1.65[%] 

600 0.98[%] 0.93[%] 1.57[%] 2.52[%] 

800 0.66[%] 0.63[%] 1.12[%] 1.84[%] 



(lineshape) with a cut p^. > 0.25 Mzz, where Mzz is the invariant mass of the ZZ-pair. 

In Figured] we show the lineshape for a Higgs mass of 600 GeV. The black line gives the full gg — > ZZ 
process at LO; the cyan line gives signal plus background (LO) neglecting interference while the blue line 
includes both gg and qq initial states (LO). The red line gives the LO signal with different cuts on the Z 
transverse momentum. 
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In Figure Owe present options for including higher-order effects. The black line is again full LO gg — > ZZ 
result, the brown line gives the multiplicative option, the red line is the additive option while the blue line 
is the intermediate option. The cyan line gives signal plus background (LO) neglecting interference. 

In Figure [3] we show the same set of results as in Figure [5] but for a Higgs boson mass of 700 GeV. 

There are different options for showing the effect of the interference, e.g. 

nNNLO J_ 

cB ~ D^(S)+D^(B) ' dM zz { ' 

However, our preferred way will be to use the following equation: 

_ KnD-°(S) + (K^) 1/2 D"°(I) + D»°(B) da 

To summarize: 

Eq. (TTfl) is our recipe for estimating the theoretical uncertainty in the effective NNLO distribution: the 
intermediate option gives the central value, while the band between the multiplicative and the additive 
options gives the uncertainty. 

Note that the difference between the intermediate option and the median of the band is always small if not 
far away from the peak where, in any case, any option becomes questionable. The ratio / K^> can be 
greater than one in some region, e.g. for MzzgtrsimSlGUGeV, almost /ih -independent with a maximum of 
1.024 at Afzz = ITeV. 

In TableH] we show the estimated theoretical uncertainty for the fractional interference correction to the 
700 GeV resonance, R' eS defined in Eg. (ITT)). The effect computed according to Eq. p"?]) is very similar to the 
one obtained by considering LO alone, shifted to the left and slightly less destructive for high values of M%%. 
In Figure H^Figure [5]) we present R' cS (Eq. ([T7| ) for = 700 GeV^(800 GeV) summarizing the percentage 
effect of interference in the effective NNLO theory. The black line gives the central value while the two blue 
lines represent the estimated theoretical uncertainty in including the NNLO if -factor. 

In Figure [5] we present the LO interference effect for = 600,700,800 GeV. In Figure [7] we present 
the effective NNLO invariant mass distribution [in — 400, 500, 600, 700, 800 GeV, including our estimate 
of the theoretical uncertainty. In Figure [8] we show the effective NNLO ZZ invariant-mass distribution for 
/in = 700 GeV including theoretical uncertainty and a comparison between 7 TeV and 8 TeV. 

In Figure [9] we present the sum signal + interference for = 400,500,600,700,800 GeV, including 
our estimate of the theoretical uncertainty. This quantity has no direct physical meaning but represents the 
pseudo-observable preferred by the experimental Collaborations. 

In Table [2] we show the effect of the pf cut; from p\ > 0.25 M zz to p\ > 0.15(0.05) M zz the signal is 
reduced by only 10%(17%) while the background is reduced by a factor 1.98(3.06). Finally in Figure [T0l we 
compare the interference effects with pj, > 0.25 M zz and with > 0.15 M z ; within the ±7h window the 
change is negligible and becomes larger for lower or higher values of A/ zz , as expected. 



Table 2: Effect of the pf 1 cu t on the total LO cross-section (Signal, Background and Total multiplied by B — 4.36 TO 
the BR for both Z bosons to decay into e or p.) for 2 Mz < Mzz < 1 TeV. 
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3.1 Residual theoretical uncertainty 



In our results we have not included uncertainties coming from QCD scale variations and from PDF +a s ; 
due to the scaling of the LO result, these uncertainties are coming from the numerator in the if -factor 
and are the typical NNLO uncertainties in gluon-gluon fusion pQ. Also excluded is the residual electroweak 
uncertainty for the signal lineshape [5] . To give an example of the complete set of theoretical uncertainties 
we select /ih = 700 GeV and define maximum and half-maxima of the signal lineshape; they are given by 

S (Mzz) = , S (Mi) = max S (M zz ) , S (mJJ = i S (Mi) . (18) 

We find Mi = 701 GeV and M~ /2 = 565 GeV, M+ 2 = 761 GeV. We define theoretical uncertainties (THU) 
according to the following sources: 1) intrinsic, the full band between multiplicative and additive options, 
2) electroweak, due to THU on 7h and Th^zz(M)zz) and described in Sect. 7 of Ref. 5 and 3) QCD 
scales. THU arise from uncertainties in underlying theoretical paradigm. Results are shown in Table [3] for 



Table 3: Theoretical uncertainty on R' eS = I/(S + B) for /ih = 700 GeV. 

~ Mz~z R' cS [%] intrinsic [%] EW[%] QCD scales 

Mi +3.82 -1.28 +1.85 -0.34 +0.26 -0.38 +0.67 

My 2 +17.19 -5.59 +8.02 -0.34 +0.44 < 0.1 

M+ /2 -10.47 -5.18 +3.58 -0.56 +4.24 -1.58 +0.82 



R' eS = I/(S + B). THU on Mi is tiny, intrinsic THU is large for the half-maxima, electroweak THU remains 
small in the window between the two half-maxima. As far as QCD scale variation is concerned we observe 
that R' eS is very stable when varying Mzz/4 < MR, A^f < Mzz- The percentage correction R' cS is the scaling 
factor that one has to apply to her/his own calculation of S + B. What it is meant in Table[3]is the following: 
the THU on I induced by QCD scale variation is between 5.39% [S + B] (M\ ss zz/4) and 6.44% [S + B] (M zz ) 
at Mzz = Mi = 701 GeV, etc The largest uncertainty in scale variation is due to the background which is 
only known at LO; at Mzz = Mi we observe a variation of 19% in Kd S + B while the variation in S + B 
is 89%. For this reason it would be difficult to work completely at LO. 



4 Conclusions 

In this paper we have addressed some issues concerning the inclusion of interference effects in gluon-gluon 
fusion, especially for high values of the Higgs boson mass. 

The results of Figure Q] - Figure [3] suggest the following compromise for effectively including higher order 
effects in the interference between Higgs and continuum contributions in gg — > ZZ. For the heavy Higgs 
scenario, above the Higgs boson peak the multiplicative (or at least the intermediate) option is recommended 
while the additive (or the intermediate) one should be preferred below the peak. However, one should also 
provide an estimate of the corresponding theoretical uncertainty. For this reason a conservative assessment of 
interference effects is represented by a central value given by the intermediate option of Eq. (|13l) with remain- 
ing theoretical uncertainty given by the full band between theadditive option of Eq.([7|) and themultiplicative 
option of Eq.flEJl. 

For an inclusive quantity the effect of the interference, with or without the NNLO if -factor for the signal, 
is almost negligible. For distributions this is radically different and we have shown our results for the ZZ 
invariant mass distribution: close to Mzz = Mh the uncertainty is small but becomes large in the rest of the 
search window [fin — 7h , fin + 7h]- The effect of the LO interference, w.r.t. LO S + B, reaches a maximum 
of +16% before the peak (e.g. at /in = 700 GeV) while our estimate of the scaled interference (always 
w.r.t. LO S + B) is 86 +7 —3 % in the same region, showing that NNLO signal effects are not negligible. 
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The estimate of the uncertainty is certainly a conservative one; however, it would be unsave to select less- 
conservative choices since we are not able to properly quantify the NLO corrections for the background. 
The percentage effect of the interference is slightly distorted when we go from LO to effective NNLO (at 
least in the intermediate option) , however the global effect on the complete distribution is sizable, as seen in 
Figure H 

In summary, we have discussed options to simulate NNLO corrections for the continuum interference 
effects to the ZZ -signal for the Standard Model Higgs boson produced via gluon fusion. The effects are large 
in the heavy Higgs-mass scenario, depending on the ZZ invariant mass. 

We can anticipate precisely what the likely criticism will be, therefore we must clearly state that LO 
kinematics is different from the NLO(NNLO) one and the if -factor will depend on aspects of kinematics 
that are not present in the LO background-interference, e.g. contribution to coming from emission of extra 
gluons etc. Therefore, there is absolutely no guarantee that different distributions will not be distorted by the 
procedure and we will need to check the effect of the NLO effects on the interference for the full kinematic 
distributions. However, one should remember that the band representing the theoretical uncertainty has 
no statistical meaning, at most a flat prior to represent maximal uncertainty. Note that the so-called 
central value also has no special meaning, although it is extremely useful for the experimental analysis where 
the MonteCarlo events are reweighted by using some analytical function. An alternative view considers 
the difference between multiplicative and additive options as a systematic uncertainty resulting in overall 
distorsion of the shape. We can try to turn our three measures of the lineshape into a continuous estimate 
in each bin; there is a technique, called "vertical morphing" [13], that introduces a "morphing" parameter / 
which is nominally zero and has some uncertainty. If we define 

D° = f a , option I, D + = maxA m D D~ = mhiA m D , (19) 
dMzz 

the simplest "vertical morphing" replaces 

D° -> D° + L (£>+ - IT) . (20) 

Of course, the whole idea depends on the choice of the distribution for /, usually Gaussian which is not 
necessarily our case; instead, one would prefer to maintain, as much as possible, the LO cancellations around 
the pealo We would like to elaborate on this with the following heuristic argument: how does the lineshape 
uncertainty band translate into uncertainty for the total cross-section? We define two curves 

C M (A , M zz ) = XDi (Af zz ) + (1 - A) D M (M zz ) , C A (X, M zz ) = A D 1 (M zz )+(1 - A) D A (M zz ) , (21) 

where Di, i = I, A, M is the lineshape according to I, A and M options. The parameter A, with < A < 1, 
parametrizes how far we are from the central value (we assume that A has a flat distribution). The uncertainty 
on the total cross-section is obtained by integrating over Mzz, once along Cm and a second time along Ca; 
the difference gives the uncertainty, maximal for A = 0. The observation is not trivial since the two curves 
cross shortly after the peak and one should not integrate over min{CA , Cm} and over max{CA , Cm}- 

It is worth noting that even for the signal alone, gg — > H — > VV, the current generation of events is done 
with something like POWHEG [H] for the initial state kinematics plus PYTHIA [45146] for the decay. The 
resulting events are rescaled in cross-section to NNLO for gg — > H while the branching ratios are rescaled 
according to Prophecy4f [24125] , For the kinematics itself this means gg — > H is NNLO while the H decay 
has the LO model of PYTHIA. 

Further study is warranted of exact NLO background for all selected channels and there is not much 
we can add at this stage; phases in NLO (NNLO) corrections to the signal are not available as well as 
NLO corrections to the background. The typical example that we have in mind is the process pp — > yyX, 
as described in Ref. [17] for the full NNLO QCD corrections; in the gg -channel the box contribution was 
computed in Ref. [48] and the next-order gluonic corrections in Ref. [49] . This is a case where all contributions 
for the gg -channel have been included in a fully-consistent manner. 

2 We gratefully acknowledge S. Bolognesi for suggesting this alternative. 
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The options that have been introduced in the literature and that we have summarized and extended 
here can only give us a rough approximation of the true result. However, this solution improves upon the 
previous estimate of 150 /Zh[%] (mh m TeV) made without any reference to the existing LO calculation of 
the interference. 

The calculation described in this work refers to a SM Higgs boson; however, ATLAS and CMS quantify 
their analysis in terms of 

Hypothesis (p,) = fia{S) + a{B) (22) 

where interference is only included in the uncertainty; /i = is the null hypothesis, fi ^= is the alternate 
hypothesis (Neyman - Pearson lemma). Another issue, currently under discussion, is the following: should 
one define the cross-section for events containing a Higgs boson in terms of sigma(S) + o~(I) with <r(B) 
defining the null hypothesis? The obvious criticism is that signal + interference is not positive defined. 
Therefore, it is not completely clear how the hypothesis should be modified to include the interference |50j . 
a possibility being [S*T] 

Hypothesis(/z) = (j, a(S) + y/fia(I) + a(B), (23) 

which is reasonably good for i.e. for analyzing the SM hypothesis. In general it would be better to 

start with an effective Lagrangian containing anomalous couplings [3] {A} (such that Aj = is the SM) and 
to define 

Hypothesis({A}) = a(S, {A}) + a(I, {A}) + a(B, {A}), (24) 
which will respect, among other things, the unitarity cancellations requested by Eq. ([¥]), i.e. 

Hypothesis(Ai) = 0, (25) 

M zz ->oo 

which proves that, in any consistent theory (not only the SM), the background knows the signal (at least 
asymptotically). Clearly, Eq. (|2"5]) imposes \i = 1 in Eq. (f2"3")l if <r(S) etc. are the SM cross-sections. A 
possibility in implementing Eq. (l2"4"|) is to use the Buchmiiller - Wyler basis |52j . 

If one does not have a clear idea of what the BSM signal is, it is difficult to optimize an analysis (Neyman 
- Pearson lemma cannot help). Alternative strategies require that, in absence of a signal, a C.L. limit is 
set on /l«h- Based on the observation that a ultra-heavy Higgs boson does not lead to a pronounced peak 
structure, Baur and Glover |18|19j made an alternative proposal for the null hypothesis: here the background 
corresponds to the minimum of S + I + B for all Higgs masses. More precisely, in ZZ -production, they define 

[N (m , mo) - N f/ijj , m J] 2 

•Smax = max5(m ) , S (mo) = -, r , (26) 

N f/4 , m J 

where Mzz > m o an d N (/^h > m o) is the number of signal events for a Higgs boson of mass /xjj. Therefore, 
Smax gives a quantitative measure of how well the /in hypothesis can be discriminated from the fi'^ hypothesis. 
The background corresponds to /i^ = but we could have /xh = 125 GeV as well. The strength of the analysis 
should be in terms of the C.L. at which the /j,h hypothesis can be excluded. Note that mo del- independent 
searches have been addressed in Ref. [53154] and recommendations have been made in Ref. [55] . 

If the signal for a light Higgs boson will be confirmed at LHC then any heavy scalar boson must be beyond- 
Standard-Modcl (BSM). As a prototype we can take a THD model; from the perspective of searching for the 
heavy partner (s) the light one is assimilable to the background (i.e. the two resonances do not interfere). Of 
course, the signal has to be rescaled according to the BSM model (couplings of the heavy partner(s), H, A, W* 1 , 
to fermions and gauge bosons) but we have not attempted a detailed analysis, e.g. taking into account the 
mass difference among the heavy bosons, related to the breaking of the custodial SU(2) symmetry. Other 
possibilities include a heavy scalar singlet with a large vacuum expectation value that can evade the potential 
instability of the SM electroweak vacuum [56 . The qualitative aspects of this work will not change. 

Finally, the interference effects for a light SM Higgs boson will follow the same pattern described in 
this work; the final state will contain four fermions (below the ZZ or WW thresholds) but the back- 
ground/interference will always be at LO and we will be missing large if -factors. Numerics will change 
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but our recipe for estimating effective NNLO signal + background + interference and the corresponding 
theoretical uncertainty will remain the same. Only progress, i.e. new more accurate calculations will be able 
to produce more accurate estimates for the theoretical uncertainty. 
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Table 4: Interference effect (percentage), R' c ^ for fin = 700 GeV with corresponding theoretical uncertainty. Here 
p| > 0.25 Mzz and 2M Z < M zz < 1 TeV. 



Bin[ GeV] 




%] 


minus error [ 


%] 


plus error [ 


%] 


210 - 212 





.24 


-0 


.07 


+0 


.11 


230 - 232 





.31 


-0 


.09 


+0 


.13 


250 - 252 





.38 


-0 


.12 


+0 


.17 


270 - 272 





.47 


-0 


.14 


+0 


.20 


290 - 292 





.56 


-0 


.17 


+0 


.24 


310 - 312 





.73 


-0 


.22 


+0 


.32 


330 - 332 





.96 


-0 


.29 


+0 


.41 


350 - 352 


1 


.47 


-0 


.44 


+0 


.63 


370 - 372 


2 


.21 


-0 


.67 


+0 


.96 


390 - 392 


2 


.82 


-0 


.87 


+1 


.25 


410 - 412 


3 


.67 


-1 


.14 


+1 


.64 


430 - 432 


4 


.72 


-1 


.48 


+2 


.13 


450 - 452 


6 


.00 


-1 


.90 


+2 


.74 


470 - 472 


7 


.49 


-2 


.37 


+3 


.41 


490 - 492 


9 


.24 


-2 


.95 


+4 


.24 


510 - 512 


11 


.30 


-3 


.63 


+5 


.21 


530 - 532 


13 


.49 


-4. 


.36 


+6 


.27 


550 - 552 


15 


.87 


-5 


.14 


+7 


.37 


570 - 572 


17 


.75 


-5 


.78 


+8 


.30 


590 - 592 


19 


.12 


-6 


.25 


+8 


.97 


610 - 612 


19 


.43 


-6 


.39 


+9 


.19 


630 - 632 


18 


.16 


-5 


.99 


+8 


.62 


650 - 652 


15 


.39 


-5 


.11 


+7 


.36 


670 - 672 


11 


.29 


-3 


.77 


+5 


.42 


690 - 692 


6 


.43 


-2 


.15 


+3 


.10 


710 — 712 


1 


no 

.28 


— 


.43 


+0 


.62 


730 - 732 


-3 


.72 


-1 


.82 


+1 


.26 


750 - 752 


-8 


.34 


-4. 


.11 


+2 


.84 


770 - 772 


-12 


.43 


-6 


.18 


+4 


.26 


790 - 792 


-16 


.05 


-8 


.06 


+5 


.55 


810-812 


-19 


.21 


-9 


.73 


+6 


.69 


830 - 832 


-21 


.98 


-11 


.19 


+7 


.69 


850 - 852 


-24 


.36 


-12 


.54 


+8 


.58 


870 - 872 


-26 


.52 


-13 


.88 


+9 


.45 


890 - 892 


-28 


.37 


-14 


.98 


+10 


.18 


910-912 


-30 


.06 


-16 


.07 


+10 


.88 


930 - 932 


-31 


.63 


-17 


.18 


+11 


.57 


950 - 952 


-33 


.01 


-18 


.10 


+12 


.16 


970 - 972 


-34 


.28 


-19 


.17 


+12 


.79 


990 - 992 


-35 


.50 


-20 


.11 


+13 


.36 
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Figure 1: The ZZ invariant mass distribution in the OFFP-scheme of Ref. [5] with running QCD scales 
for (Liu = 600 GeV. B — 4.36 ■ 10 -3 represents the BR for both Z bosons to decay into e or |i. The black 
line gives the full gg — > LL process at LO; the cyan line gives signal plus background (LO) neglecting 
interference while the blue line includes both gg — > ZZ and qq — > TIL components (LO). The red line 
gives the LO signal. 
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Figure 2: The ZZ invariant mass distribution in the OFFP-scheme of Ref. [5] with running QCD scales 
for fiu = 600 GeV. B — 4.36 ■ 10~ 3 represents the BR for both Z bosons to decay into e or (i. The black 
line is the full LO gg — > TIL result, the brown line gives the multiplicative option of Eq.([8|, the red line 
is the additive option of Eq.0 while the blue line is the intermediate option of Eq. (ll3[) . The cyan line 
gives signal plus background (LO) neglecting interference. 
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Figure 3: The ZZ invariant mass distribution in the OFFP-scheme of Ref. [5] with running QCD scales 
for fiu = 700 GeV. B — 4.36 ■ 10~ 3 represents the BR for both Z bosons to decay into e or (i. The black 
line is the full LO gg — > TIL result, the brown line gives the multiplicative option of Eq.([8|, the red line 
is the additive option of Eq.0 while the blue line is the intermediate option of Eq. (ll3[) . The cyan line 
gives signal plus background (LO) neglecting interference. 
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Figure 4: Interference effects (see Eq. (|17[0 in the ZZ distribution for fm 
the central value, the blue lines give the estimated theoretical uncertainty. 



700 GeV. The black line is 
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Figure 5: Interference effects (see Eq. (|17[0 in the ZZ distribution for [in — 800 GeV . The black line is 
the central value, the blue lines give the estimated theoretical uncertainty. 
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Figure 6: LO Interference effects [%] in the ZZ distribution for /m = 600, 700, 800 GeV. 
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Figure 7: Effective NNLO ZZ invariant-mass distribution for fin = 400, 500, 600, 700, 800 GeV including 
theoretical uncertainty. B — 4.36 • 10 -3 represents the BR for both Z bosons to decay into e or |i. 
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Figure 8: Effective NNLO ZZ invariant-mass distribution for fin = 700 GeV including theoretical 
uncertainty and a comparison between 7 TeV and 8 TeV. B = 4.36 ■ 10 -3 represents the BR for both Z 
bosons to decay into e or p.. 
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Figure 9: Effective NNLO ZZ invariant-mass distribution for fi H = 400, 500, 600, 700, 800 GeV including 
theoretical uncertainty. Only signal + interference is plotted. B = 4.36 • 10 -3 represents the BR for 
both Z bosons to decay into e or p.. 
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,mre 10: Interference effects (see Eq.JT7|)) in the ZZ distribution for /ih = 700 GeV comparing 
> 0.25 M zz (blue) with p| > 0.15 M zz (red). 
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